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by 


The  equations  of  motion  of  a deformable  aircraft  are  developed  in  detail 
from  Lagrange's  equations  for  a non-inertial  frame.  Particular  account  is  taken 
of  the  influence  of  the  propulsive  and  effective  forces  produced  by  power  units 
containing  rotating  parts.  The  development  ventures  to  a certain  extent  into 
the  non-linear  regime.  By  an  appropriate  choice  of  deformation  modes  the 


principal  frame  of  reference  can  be  specified  as,  for  example,  mean-body  axes  or 
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1 INTRODUCTION 

There  have  been  a number  of  recent  papers  on  the  development  of  equations  of 
motion  of  deformable  aircraft  ( eg  Refs  1-5  and  8-10).  One  may  well  ask  isn't 
that  enough?  To  help  the  reader  form  his  opinion,  we  state  here  the  purpose  of 
the  present  paper.  It  is  briefly,  to 

(i)  derive  in  detail  the  equations  of  motion  of  a deformable  aircraft 
using  an  arbitrary  non-inertial  frame  as  the  frame  of  reference  for 
Lagrange's  equations*; 

(ii)  develop  these  equations  up  to  and  including  all  second  order  terms 
in  the  deformational  and  body  freedom  coordinates; 

(iii)  take  account  in  the  equations  of  motion  of  the  fact  that  part,  or 
parts  of  the  propulsive  unit  may  be  rotating. 

In  the  past  when  a non-inertial  frame  of  reference  has  been  used,  it  has 
almost  without  exception  been  either  body-fixed  axes  (origin  and  orientation 
fixed  in  a small  portion  of  the  aircraft)  or  mean-body  axes  (origin  at  the  centre 
of  gravity  and  orientated  so  that  the  kinetic  energy  relative  to  the  axes  is 
a minimum) . The  deformation  modes  had  therefore  to  satisfy  certain  particular 
conditions.  In  Ref  2,  the  equations  of  motion  were  in  effect  obtained  for  an 
arbitrary  non-inertial  frame  by  transformation  from  those  obtained  using  body- 
fixed  axes  but  they  are  slightly  different**  from  those  here  obtained  as  a conse- 
quence of  a different  representation  of  the  deformations.  Here  we  assume  the 
deformations  relative  to  the  arbitrary  non-inertial  frame  are  precisely  a linear 
function  of  the  generalised  coordinates  while  in  Ref  2 it  was  just  the  deforma- 
tions relative  to  the  body-fixed  axes  which  had  this  characteristic.  In  Ref  2 
the  forms  of  perturbation,  when  using  either  of  the  non-inertial  frames, 
arbitrary  or  body-fixed,  were  selected  so  that  any  perturbation  which  could  be 
described  by  the  one,  could  also  be  precisely  described  by  the  other. 

Non-linear  effects  is  a convenient  phrase  when  one  has  discrepancies 
between  predictions  and  experience,  but  few  have  made  any  attempt  to  include 
non-linearities  in  their  mathematical  model,  and  when  they  have,  it  has  only 
been  in  a restricted  way  - non-linear  behaviour  of  a power  control  unit,  for 
example.  To  aid  any  more  general  non-linear  investigation  it  did  therefore  seem 

* Alternatively  one  can  think  of  this  as  'using  a non-inertial  frame  ...  in  con- 
junction with  an  arbitrary  set  of  deformation  modes. 

**  The  linear  approximations  to  tha  equations  of  motion,  though  not  the  same,  are 
exactly  equivalent  since  the  linear  approximations  to  the  perturbations  are 
the  same. 
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worthwhile  to  carry  out  the  present  development  of  the  equations  of  motion  as  far 
as  the  second  order  terms  in  the  generalised  coordinates.  However,  before  one 
performs  any  non-linear  dynamical  analysis  one  should  be  firmly  convinced  that 
its  results  will  be  meaningful.  Solution  will  not  be  easy  in  general;  and  one 
should  ask  oneself  "Is  my  data  adequate?  Are  my  assumptions  about  the  form  of 
the  local  forces  sufficiently  near  the  truth?"  and  so  on  ( of  Ref  4).  In 
particular,  in  connection  with  this  work  one  may  well  ask  what  possibly  can  be 
gained  by  the  use  of  the  second  order  equations.  They  certainly  won't  make  any 
change  in  the  stability  or  otherwise  of  the  datum  state.  They  may,  however, 
predict  a boundary  to  the  extent  of  the  stability,  whereas  with  the  linear  system 
if  we  have  stability  then  it  is  always  stability  'in  the  large'  (and  vice  versa). 
There  is  little  value  in  having  a stable  datum  state  if  the  extent  of  the 
stability  is  very  small,  and  conversely  an  'unstable'  datum  state  may  be  of  little 
concern  if  the  extent  of  the  instability  is  very  small. 

Years  ago  someone  made  the  assumption  "Oh,  we  can  neglect  that"  - perhaps 
because  it  didn't  seem  worth  the  effort  to  do  anything  else  - and  their  successors 
have  followed  in  their  footsteps  without  question,  or  so  it  appears.  One  such 
assumption  was  to  take  no  account,  in  aircraft  dynamical  studies  such  as  stability 
and  flutter,  of  the  fact  that  the  aircraft  propulsive  unit  contains  one  or  more 
rotating  parts  having  an  appreciable  moment  of  inertia.  Consideration  is  some- 
times given  to  non-linear  aerodynamic  terms  which  couple  symmetric  and  anti- 
symmetric freedoms;  but  may  not  the  linear  gyroscopic  coupling  between  the  same 
freedoms,  resulting  from  engine  rotation,  be  just  as  important?  And  there  lies  the 
the  reason  for  the  third  aspect  of  the  present  development. 

These  last  two  aspects  - non-linearities  and  rotating  parts  - add  con- 
siderable complication  to  the  development  (of  eg  the  size  of  section  7).  Thus, 
partly  for  the  benefit  of  those  who  are  mainly  interested  in  the  first  purpose  of 
the  present  paper,  and  also  to  show  how  other  developments  of  the  equations  of 

motion  can  be  considered  as  particularisation  of  the  development  using  an 

. . 12  ... 
arbitrary  non-inertial  frame,  a second  paper  has  been  written  omitting  the 

aspects  (ii)  and  (iii). 

2 THE  SEMI-RIGID  MODEL 

We  wish  to  derive  equations  of  motion  which,  to  a good  approximation, 
govern  perturbations  of  the  motion  of  an  aircraft,  from  a datum  state  in  which 
it  is  flying  with  constant  linear  velocity  and  zero  angular  velocity  in  a uniform 
atmosphere.  As  a basic  frame  of  reference  we  will  take  a 'constant-velocity' 
frame  which  is  such  that,  during  the  datum  motion,  the  aircraft  is  stationary 
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relative  to  it.  In  addition  a second  frame  of  reference  is  used  which  is  such 
that  the  aircraft  may  be  considered  as  fixed  in  this  frame  when  undeformed,  even 
during  the  perturbed  motion.  This  frame  is  identified  therefore  by  a set  of  axes 
called  the  ' no-deformation-body-fixed * axes';  and  deformations  are  considered, 
by  definition,  to  be  departures  from  the  positions  defined  to  be  the  undeformed 
positions  within  this  reference  frame.  It  will  be  seen  therefore  that  the  no- 
deformation-body-fixed axes  can  in  effect  be  arbitrarily  chosen  by  the  analyst. 

In  other  words  in  going  from  the  unperturbed  to  the  perturbed  position  of  the 
aircraft  he  puts  down  a stepping  stone  inbetween  just  where  he  likes  - she  (the 
aircraft)  steps  rigidly  to  the  stone  and  from  there  she  deforms  herself  to  make 
the  second  step  to  the  perturbed  position.  It  should  be  appreciated  however  that 
the  choice  of  a set  of  deformation  modes  does,  in  effect,  select,  in  the  perturbed 
state,  an  'optimum'  position  for  the  no-deformation-body-fixed  axes.  Thus,  for 
example,  these  axes  become  mean-body  axes  if  the  modes  satisfy  certain  conditions*, 
and  they  become  body-fixed  axes  when  the  modes  satisfy  another  set  of  conditions*. 

Let  the  constant-velocity  frame  of  reference  be  identified  by  constant- 

velocity  axes  0 x y zc  which  are  such  that,  in  the  datum  motion,  they  coincide 

with  some  axes  fixed  in  the  aircraft.  The  location  of  the  no-deformation-body- 

fixed  axes  (0  x y z ) can  then  be  specified  by  the  translations  and  rotations 
n n n n 

which,  when  applied  to  the  constant-velocity  axes,  would  make  the  two  sets  of 
axes  coincident.  Let  these  translations  and  rotations  be,  in  succession: 

(i)  Translations  XjC\  yjC\  i-n  the  directions  of  the  respective 

constant-velocity  axes. 

(ii)  Successive  rotations  ij>,  8,  $ about  the  carried  axes  °nz*  0ny*  ®nx  • 

These  translations  and  rotations  will  be  taken  as  the  six  'body-freedom' 
generalised  coordinates.  They  uniquely  specify  the  position  of  the  no- 
deformation-body-fixed axes. 


* These  two  sets  of  conditions  are, 
later  (equation  ( 2-1 )). 

(i)  ^ 6mR  • V 6mA  R - 0 


L-j 


l—i 


in  terms  of  the  modal  matrix 
for  mean-body  axes. 


R introduced 


(ii)  At  Xj  • y^  ■ ■ 0 , R ■ 0 , has  null  second  and  third 

dR  dR 

rows,  t — has  null  first  and  third  rows,  -r—  has  null  first  and 

*yf 

second  rows  - for  body-fixed  axes. 
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We  will  assume  Chat  the  second  step,  the  deformations,  can  be  written  as  a 
linear  combination  of  n deformational  degrees  of  freedom;  and  so  the  position 
of  a particle  relative  to  the  no-deformation-body-fixed  axes  is  given  by 


r (to 

X 

n 

- 

xf 

+ R 

ql 

(n) 

yn 

yf 

e 

• 

e 

Z(n) 

si 

ln 

_zf 

qn 

(2-1) 


where  R is  a modal  matrix  whose  elements  are  functions  of  the  unperturbed 
coordinates  (x^,  y^,  z^)  of  the  particle.  Note  that,  in  doing  so,  we  are  not 
making  any  assumption  that  the  deformations  are  small.  We  are  however  excluding 
the  exact  representation,  in  one  degree  of  freedom,  of  the  rotation  of  a part  of 
the  aircraft  as  a rigid  body,  but  this  should  be  of  minor  importance  ( cf  the 
treatment  of  engine  rotation  in  section  7). 


Thus  the  position  of  a typical  particle  relative  to  the  origin  of  the 
constant-velocity  axes  and  referred  to  those  axes  is  ( of  Appendix  A of  Ref  1) 
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and  is  such  that 
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If  now  we  substitute  in  (2-2)  the  expansion  of  S 
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where  J,  is  defined  by  equation  (3-4). 
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3 VELOCITIES 

Writing  the  linear  velocity  components  of  the  constant-velocity  axes, 
referred  to  those  axes,  as  |u^,  v^,  w^f,  and  remembering  that  their  angular 
velocity  is  zero,  it  is  easily  seen  ( af  Ref  1)  that  the  linear  and  angular 
velocities  of  the  no-deformation-body-fixed  axes  are 
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Consequently,  using  (2-8), 


the  linear  velocity  of  a particle  is  given  by 
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The  above  expression  assumes  that  the  modal  matrix  R and  the  unperturbed 
coordinates  (x^,  y^,  z^)  are  constant.  However  when  we  subsequently  come  to  con- 
sider the  effect  of  a rotating  part  of  an  engine  (section  7.1)  this  will  not  be 
so.  The  appropriate  expression  for  the  velocity  of  a particle  then  is 
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This  expression  will  be  required  subsequently  (section  4.2)  in  determining  the 
effective  forces. 
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4 EQUATIONS  OF  MOTION 

4. 1 For  the  deformational  degrees  of  freedom 

The  non-inertial  form  of  Lagrange's  equation  is*'^ 


3V0 

— - + J.  + G. 
3q.  l l 
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•Hi) 


3W 

3«i 


" Qi 


(4-1) 


where  V 
w 
G. 
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J. 
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is  the  centrifugal  potential  function 

is  the  kinetic  energy  relative  to  the  frame  of  reference 
is  the  gyrostatic  force 

is  a rather  similar  force  to  G^  but  depending  on  the  angular 
accelerations  rather  than  the  angular  velocities  of  the  frame  of 
reference 

is  the  generalised  force  obtained  by  the  method  of  virtual  work  - in 
the  assumed  infinitesimal  virtual  displacement  the  frame  of  reference 
is  regarded  as  stationary 

q.  is  the  generalised  coordinate  of  the  ith  degree  of  freedom  where  that 
freedom  is  one  such  that  the  position  and  orientation  of  the  frame  of 
reference  is  independent  of  it. 

Thus,  this  equation  can  be  used  for  the  degrees  of  freedom  i = 1,  . ..,  n when 
taking  the  frame  of  reference  to  be  the  no-deformation-body-fixed  axes.  The 
various  quantities  listed  above  are  given  by  the  following  expressions,  where  6m 
is  the  mass  of  a particle  of  the  aircraft,  and  the  summations  are  for  all 
particles : 
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Using  the  expressions  for  the  velocities  of  no-deformation-body-fixed  axes 
(equations  (3-1)  and  (3-2))  and  for  the  location  of  a particle  relative  to  these 
axes  (equation  (2-1))  we  find  that  the  centrifugal  potential  function  is  given, 
after  some  manipulation  (of  Ref  4,  Appendix),  by 
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where  I 


»(-z  6mA^  ^ is  the  matrix  of  moments  and  products  of  inertia  of  the 


Lu  "f 

aircraft  (in  the  datum  motion  state). 
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From  (2-1),  (3-2)  and  (3-3)  the  gyrostatic  forces  are  easily  seen  to  be 


!^6mRTA«R  qj  + ... 


(4-11) 


Similarly,  the  forces  (equation  (4-5))  are 
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The  kinetic  energy  relative  to  the  frame  of  reference  (the  no-deformation-body- 
fixed  axes)  is  from  (2-1)  and  (4-3) 


and  so. 
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Consequently,  collecting  together  the  various  contributions,  we  find  that  the 
column  vector  of  the  effective  forces  (the  left-hand  side  of  equation  (4-1))  for 
the  degrees  of  freedom  1 -*■  n , is 
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Expressions  for  the  various  contributions  to  the  generalised  forces  (the  right- 
hand  sides  of  equation  (4-1))  are  derived  below  (sections  5.1,  6.1,  7.1  and  8.1). 

4.2  For  the  body- freedoms 

To  obtain  the  additional  equations  required  to  form  a complete  set  we  make 
use,  as  in  Refs  1 and  2,  of  the  principles  of  linear  and  angular  momentum  applied 
respectively  in  the  directions  of  and  about*  the  no-deforraation-body-f ixed  axes. 
These  give**  (see  equation  (3-7)) 


* Strictly  speaking  about  stationary  axes  which  are  instantaneously  coincident 
with  the  no-deformation-body-fixed  axes. 


**  We  have  used  the  relationship  2B 
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Expressions  for  the  various  contributions  to  the  overall  forces  and  moments  are 
derived  below  (sections  5.2,  6.2,  7.2  and  8.2). 
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5 AERODYNAMIC  FORCES 

The  boundary  condition  to  be  satisfied  by  the  air  motion  at  the  surface  of 
the  body  involves  the  velocity  of  the  surface  normal  to  itself.  The  velocity 
vector  of  a particle  is  given  by  equation  (3-5),  when  referred  to  the  no- 
defonnation-body-fixed  axes;  and  the  slope  vector  of  the  body  surface,  referred  to 
the  same  axis,  will  ( of  equation  (2-1))  be  merely  a function  of  the  generalised 
coordinates  qj  -*■  q . The  local  aerodynamic  force  vector  will  therefore  be  a 
function  of  the  three  quantities 


(5-1) 


■ 

• (cf 

• 

- -| 

uf 

+ 

xi 

• % 

& 

9 

ql 

• (c) 

• 

a 

vf 

yl 

► 

0 

• 

• 

• (c) 

9 

k 

w 

Ld 

zl 

4> 

s 

1 —4  1 1 

and  their  derivatives  with  respect  to  time.  As  in  previous  papers  ’ we  will 
take  a model  containing  no  hereditary  constituent.  Even  so,  it  is  difficult  to 
decide  what  second  order  approximation  one  should  take  for  the  local  aerodynamic 
force  vector.  The  linear  approximation  should  clearly  have  the  form,  bearing  in 
mind  (2-8)  and  (3-3), 
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(5-2) 


where  the  coefficients  e.,  e^,  e^  etc,  may  be  differential  operators.  The 
problem  is  in  the  choice  of  the  second  order  terms.  If  there  were  no  deforma- 
tions one  could,  with  a little  superfluity,  cater  fairly  adequately  for  such 
terms  by  the  expression  ( of  equations  (2-9),  (3-4),  (4-7)  and  (4-8)) 
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(5-3) 


where* 


E(uv)  _ E(vw) 


(5-4) 


(5-5) 


and  the  elements  of  the  E matrices  may  be  differential  operators.  We  note  that, 
if  P and  Q are  constant  3 * 3 matrices,  then 
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In  addition  to  (5-3),  we  can  probably  cover  any  second  order  terms  involving 
products  of,  on  the  one  hand,  a body  freedom  coordinate,  and  on  the  other  a 
deformational  coordinate,  by  an  expression  of  the  form 


* A term  E^' 
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where  the  E^  matrices  are  3><3  matrices  whose  elements  may  be  differential 
operators,  and  the  matrices  are  constant,  3 x n,  matrices.  Further  terms 

on  similar  lines  could  be  added  to  cater  for  any  second  order  terms  involving  the 
product  of  two  deformational  coordinates.  We  will  certainly  not  do  this,  however. 
To  do  so  would  necessitate  the  inclusion,  for  consistency,  of  similar  terms  in 
the  expression  for  the  local  structural  force  (see  section  8);  and  the  two 
together  would  require  a mass  of  coefficients  whose  prediction  or  measurement 
seems  a faint  possibility  in  the  dim  and  distant  future.  As  it  is  (5-3)  and 
(5-7)  introduce  a further  (351  + 45n)  matrix  elements  beyond  the  (21  + 3n) 
elements  of  the  linear  approximation.  The  contributions  from  these  terms  to  the 
overall  and  generalised  forces  moreover  admit  of  no  simplification  after  their 
original  specification  (of  Ref  4,  section  4).  We  will  therefore  in  the  two  subse- 
quent sections  (5.1  and  5.2)  derive  the  forces  resulting  from  the  approximation 
to  the  local  aerodynamic  force  vector  which  is  linear  in  the  three  quantities 
(5-1),  and  leave  any  extension  to  the  reader  if  he  so  desires.  This  approximation 
is 
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where  we  have  made  use  of  the  relationship 
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With  the  expression  (5-8)  for  the  local  aerodynamic  force  vector,  the 
column  vector  of  the  generalised  aerodynamic  forces  is 


The  translational  force  on  the  aircraft,  and  the  moment  about  the  origin 
of  the  no-deformation-body-fixed  axes,  resulting  from  the  aerodynamic  load 
distribution  (5-8),  are  respectively 
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where 
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and 
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GRAVITATIONAL  FORCES 


On  each  particle  of  the  aircraft  there  will  be,  according  to  the  standard 

assumption,  a gravitational  force  6mg  acting  in  the  z direction  (vertically 

7 . ^ 

downwards)  of  the  normal  earth-fixed  axes  . If  the  orientation  transformation 
matrix  which  takes  us  from  these  axes  to  the  constant-velocity  axes  is  S^,  , 

where 


s*£  * 


(6-1) 


(of  equations  (2-4),  (2-5)  and  (2-6));  then  the  matrix  which  performs  the 
orientation  transformation  from  the  normal  earth-fixed  axes  to  the  no-deformation- 
body-fixed  axes  is  . Thus  the  local  gravitational  force  vector  referred  to 

the  no-deformation-body-f ixed  axes  is 
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where 


- sin  0^ 
sin  cos 

cos  cos  8^ 


(6-3) 


is  the  last  column  of  . The  angles  0^,  that  we  have  introduced  in 

this  section  are  respectively,  the  angle  of  bank,  the  angle  of  inclination,  and 
the  nose-azimuth  angle,  in  the  datum  motion.  Since  our  aatum  motion  is  flight 
in  a straight  line,  there  is  no  loss  in  generality  in  taking  ^ to  be  zero; 


and  also,  for  conventional  aircraft  under  normal  flying  conditions,  the  angle  of 
bank  will  also  be  zero.  If  the  velocity  of  the  aircraft  referred  to  normal 

earth-fixed  axes  is  { U£Q » vfQ»  wfo^’  during  the  datum  motion,  then  its  velocity 
referred  to  the  constant-velocity  axes  is 


uf  = ufo4f  + vfOk*f  + wf0*4>f 
vf 


(6-4) 


where,  for  ^ » 0 , 
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6. 1 Generalised  gravitational  forces 

The  column  vector  of  the  generalised  gravitational  forces  is 
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6.2  Overall  gravitational  forces 

Using  the  expression  (6-2)  for  the  local  gravitational  force  vector,  it  is 
easily  seen  that  the  overall  translational  forces,  and  moments  about  the  origin 
of  the  no-deformation-body-f ixed  axes,  are,  referred  to  those  axes,  given  by 


W J 

(6-10) 


7 PROPULSIVE  FORCES 

In  Refs  2,  3 and  4,  a very  simple  model  of  the  propulsive  forces  was  used; 
and,  in  particular,  no  account  was  taken  of  the  fact  that  most  aircraft  propulsive 
units  contain  one  or  more  rotating  parts  which,  in  effect,  exert  significant 
'forces*  on  the  system.  The  opportunity  is  therefore  being  taken  in  the  present 
development  to  include  a more  general,  more  sophisticated,  and,  it  is  hoped, 
potentially  more  accurate  model.  The  development  will  be  written  assuming  the 
engine  consists  of  just  one  rotating  part  but  the  extension  to  one  with  a number 
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of  rotating  parts,  or,  indeed,  to  a number  of  engines,  is  obvious  and  merely 
laborious . 

7. 1 The  model 

We  will  assume  that  the  rotating  part  is  axisymmetric  and  almost  rigid,  and 
that  in  the  unperturbed  state  it  is  rotating  with  angular  velocity  p^  about  its 
axis  of  symmetry. 

Let  two  points  on  the  axis  of  symmetry  be 
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such  that  the  direction  of  rotation  is  right-handed,  about  the  line  from  point  1 

to  point  2,  if  p is  positive;  and,  with  x , = x,»  - xr,  etc,  let  us  write 

p pt  tz  rl 


(7-2) 
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Thus  (xp,  yp,  Zp)  specify  the  orientation  of  the  axis  of  the  rotating  part.  We 
will  assume  that 
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(7-3) 


Since  we  are  assuming  the  rotating  part  is  almost  rigid  - to  make  it  completely 
rigid  would,  amongst  other  things,  require  Rj  “ Rj  which  is  an  unreasonable 


limitation*  - we  make 
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* This  is  a consequence  of  our  use  of  a strictly  linear  expression  for  the  defor- 
mations (equation  (7-1)  above  for  example)  which,  while  giving  appreciable 
simplification  does  exclude  the  exact  representation  of  a rotation  of  a part 
of  the  aircraft  as  a rigid  body. 
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Then 

(xpf  + ypf  + Zpf)  + [ql  qJ (R2  • R1)T(R2 

(7-5) 

and  so  che  perturbation  in  the  length  of  the  axis  of  the  rotating  part  between  our 
two  reference  points  is  of  second  order  in  the  generalised  coordinates. 

The  axes  transformation  matrix  (of  equations  (2-3)  to  (2-6))  which  achieves 
the  orientation  transformation  from  the  no-deformation-body-fixed  axes  to  axes 
fixed  (almost)  in  the  rotating  part  with  the  x-axis  parallel  to  the  axis  of 
synmetry  (axis  of  rotation)  will  then  be  Rp^FI  where  from  the  form  of  Rppt 
(equation  (2-4)),  II  must  be  such  that 
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where,  here  and  subsequently,  we  take 


to  mean  the  positive  value  of  the 


square  root.  Consequently,  since  we  also  require  n 
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(of  equation  (2-7)), 
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* There  is  an  arbitrariness  about  the  signs  of  the  last  two  terms  in  the 

diagonal  matrix.  The  one  chosen  ensures  that,  when  zp  ■ 0 , the  II  transfor- 
mation represents  just  one  rotation.  When  yp  ■ 0 it  can  also  represent  just 
one  rotation  if  xp  is  positive.  Otherwise  it  represents  two  rotations  and  no 

2 2 . . 

more.  One  could  use  (7-7)  when  x + y *0  but  one  then  has  arbitrarily  to 
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Expanding  n as  a Taylor  series  about  its  unperturbed  value  we  find  that, 
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where  the  B ^ are  symmetric  matrices  such  that  the 
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Expressions  for  the  various  derivatives  of  are  obtained  in  the  Appendix. 

Equation  (7-8)  can  be  written  in  terms  of  the  generalised  coordinates  by  the  use 
of  (7-2). 

The  position  of  a general  point  on  the  rotating  part,  relative  to  the 
point  a * 1 of  (7-1),  and  referred  to  the  axes  'fixed'  in  the  rotating  part, 
will  therefore  have  the  form 
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Our  assumption  of  almost  rigidity  requires  that  this  expression  contains  no  first 
order  terms  in  the  generalised  coordinates.  Let  us  assume  that  there  is  a matrix 
B relating  the  modal  function  for  a general  point  on  the  rotating  part  with  that 
for  a point  on  its  axis  of  rotation  by  the  equation 
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which  will  achieve  this.  Substituting  (7-8)  and  (7-10)  in  (7-9),  making  use  of 
(7-2),  and  equating  the  first  order  terms  to  zero,  we  obtain 
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where 
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are  the  unperturbed  coordinates  of  a particle  on  the  rotating  part,  relative  to 
the  point  a ■ 1 of  (7-1),  and  referred  to  the  no-deformation-body-fixed  axes. 
We  also  of  course  require  B to  be  such  that,  at  the  point  a = 2 of  (7-1), 
equation  (7-10)  gives  R = • Now  from  (7-12)  and  the  results  obtained  in  the 

Appendix  we  find  that  B at  this  point  is  not  the  unit  matrix  but  (of  also 
equation  (7-16),  (7-50),  (7-50a)  and  (7-50b)) 
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Thus,  the  condition  (7-4),  already  stated,  is  also  necessary;  and  so,  if  this  is 
satisfied,  the  desired  almost  rigidity  is  achieved  by  requiring  that  the  modal 
matrix  on  the  rotating  part  satisfies  (7-10)  where  B is  given  by  (7-12).  It  is 


shown  in  the  Appendix  that  each  of  the  products 
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a skew  syranetric  (A  type)  matrix.  The  same  is  therefore  true  of  the  products 
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-z n etc,  and  so  (7-12)  can  be  rewritten  as 
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Furthermore,  making  use  of  the  properties  of  A type  matrices  and  writing  (of  the 
Appendix) 
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we  have 
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The  position  of  a general  point  on  the  rotating  part,  relative  to  the 
origin  of  the  no-deformation-body-fixed  axes  and  referred  to  those  axes  thus, 
from  (7-10)  and  (7-13),  has  the  form 
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Note  that  the  x^  etc,  and  hence  B , are  not  constant  with  respect  to  time 
for  ( of  (7-9)) 
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Also  from  equation  (7-17) 
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and  so  using  these  equations  along  with  (7-17),  (7-19),  (3-1)  and  (3-2),  we  find 
that  the  contributions  to  the  effective  forces,  the  left-hand  side  of  Lagrange's 
equation  (equation  (4-1)),  from  the  rotating  part  of  the  engine,  are  given  by 
( of  equations  (4-2)  to  (4-5)) 
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Now  it  is  clear  from  the  symmetry  of  the  rotating  part  that  none  of  the  above 

functions  (equations  (7-27)  to  (7-31))  can  depend  on  the  specific  value  of  p^t 

it  doesn’t  matter  how  far  round  the  part  has  got  in  its  rotation;  and  so  the 

value  of  any  term  which  does  not  involve  a derivative  of  B or  ix^*P  z^H 

l pf  Pf  Pf  / 

will  be  independent  of  p^  . This  can  be  demonstrated  for  example  for  such 
terms  as 


StaBA^B  etc 


(7-32) 


where  n ?}  is  any  vector  which  is  constant,  with  respect  to  location,  over  the 


It  follows  that 


>(p)  . 
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(7-33) 


6m (SST  + BBT)  = 0 


(7-34) 


Thus,  using  the  relationships  (3-2)  and  (7-23) 


c -N 


f’n]  LGnJ  Pp=0 


= 2(R2  - R,)  x 


• (n) 

SmBA  , x xu‘ 
(n)  pf 

’ iff 


,./) 


(R2  ■ R])  9, 


(7-35) 


and  so  the  additional  contribution  to  the  effective  forces  due  to  the  fact  that 
part  of  the  engine  is  rotating  is 
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since 
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P P 


“ P I diagjo  1 1} 


(7-42) 


and  so  from  (7-22) 
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and  so,  using  (7-22)  and  (7-21)  along  with  these  last  two  equations 
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This  relationship  can  be  alternatively  written  as  (see  the  Appendix,  equal 
(A-23)) 
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Now  it  follows  from  equations  (7-4),  (7-10)  and  (7-16)  that* 

■a 


cntA  = i - 
f x , 
pf 


[v  ypf  Vi 


* We  do  not  require  the  vector  a 6 However,  it  can  be  shown  that 


- cn  a 

f X 


- 1 


pf 


l 


T x 

<<■  * y;„  * <,)  pf 


2 2 2 

;pf  + yPf  + J 


and  so 


a 

e 

[xlt  * ypf  * v) 

Xpf 

ypf 

Y 

LvJ 

ion 


(7-47) 


(7-48) 


1 

| 

(7-49)  I 


(7-50) 


(7-50a) 


(7-50b) 


44 


and  so  (7-49)  can  be  rewritten 
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Furthermore,  since 
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and 
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we  have,  from  (7-16)  and  (7-19) 


BB*  = - C|RT  A , .R  + RT  A , .R  |(T 

p t x(p)  p t P t x(p)  p t 

P pf  P p pf  p 


T I • *T 

- CRl  I A / vR  R A , + 

P tl  ^(p)  P„c  P„t  „<P) 

p \ xpf  p p Xpf 


A , .R  R 
X<P)  P t P t 

. pf  P P 


•T  .2 


+ (R  R A , x > + R R A , XR  R R C* 
p t P t (p)  P t p t V(P)  p t p t I pt 

P p xpf  J p p Xpf  p p / p 


2 T 
- p CR 

P Ppt 


_ Jp)  _„(P),(P) 


V 


7pf  Zpf 


-JPhW 

Pf 


ypf  V 


(7-54) 


and  so  (af  equations  (7-21)  and  (7-22)) 
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Now  from  (7-48),  differentiating 
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and  it  can  be  shown  that  (see  (7-7)  and  (7-1^ 
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Thus,  the  last  five  equations  combine 
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Z(P)  ,T 

6mBA  , required  in  equation  (7-39),  is 

perhaps  the  most  troublesome. 

In  the  datum  motion  the  velocity  of  a particle  on  the  rotating  part  relative 
to  the  axis  of  rotation  and  referred  to  the  axes  'fixed*  in  the  rotating  part,  is 
by  the  usual  formula 


-z(p) 

Pf 


(7-61) 


since  the  x^p  are  constant.  Consequently,  transforming  to  the  no-deformation- 
body-fixed  axes,  we  have 
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Using  this  equation,  along  with  (7-18)  and  (7-16),  we  have 
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Now,  as  mentioned  earlier,  the  sum  of  this  term  (*  6m)  over  the  rotating  part 
must  be  independent  of  the  value  of  p^t  , and  so 
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(7-64) 
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(7-66) 

Thus,  taking  account  of  (7-21)  and  (7-22),  as  well  as  (7-65), 
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The  additional  contribution  to  the  effective  forces  due  to  the  fact  that  part 
of  the  engine  is  rotating  is  therefore,  from  (7-39),  using  the  results  (7-43)  and 
either  (7-46),  (7-68),  (7-55),  (7-56)  and  (7-49),  or  (7-47),  (7-70),  (7-60), 
(7-51)  and  (3-2),  accordingly  as  we  wish  to  express  the  results  in  terms  of  the 
matrix  C or  not. 
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The  latter  form  has  only  been  given  explicitly  as  far  as  the  first  order  terms 
in  the  generalised  coordinates  (in  particular  the  contribution  from  (7-70)  does 
not  appear),  for  it  is  thought  the  former  form  will  be  the  more  convenient  if  one 
wishes  to  venture  further.  Thus  the  second  order  terms  in  (7-71)  are  ( af  (3-2), 
(3-3)  and  (7-88)  or  (A-23)) 
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7.3  Contributions  to  the  effective  forces  in  the  body-freedom  equations 

19 

The  right  hand  sides  of  equations  (4-17)  and  (4-18)  are  respectively  ’ 
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where  the  constituents  are  given  by  equations  (3-6),  (3-2)  and  (2-1).  Now  on 
the  rotating  part  (2-1)  becomes  (7-17)  and  its  derivative  is  (7-19).  Thus  we 
have,  from  (3-6), 
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Now,  making  use  of  the  axisymnetric  properties  (equations  (7-20)  to  (7-22))  of 
the  rotating  part  it  can  be  shown,  after  some  manipulation,  that 
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Consequently,  using  (3-2)  and  (2-13), 
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In  Ref  2 it  was  assumed  that  the  propulsive  force  acting  on  any  particular 
particle  had  constant  components  in  the  direction  of  the  body-fixed  axes*. 
However,  as  indicated  in  Ref  4,  it  is  only  convenient  to  do  this  when  the  body- 
fixed  axes  are  constrained  to  remain  mutually  at  right  angles  during  small  per- 
turbations. Body-fixed  axes  do  not  feature  at  all  in  the  present  development 
in  any  other  respect  and  so  it  seems  rather  artificial  to  introduce  them  in  a 
constrained  form  now.  We  have  however  (of  section  7.1)  used  axes  which  are 
virtually  fixed  in  the  rotating  part  (or  parts)  of  the  propulsive  unit  and  these 

axes  are  such  that  they  are  always  mutually  perpendicular.  It  is  therefore  as 

2 

good  an  assumption  as  that  used  previously  to  say  that  the  local  propulsive 
force  on  the  rotating  part  has  constant  components  in  the  directions  of  these 
latter  axes,  while  on  the  adjacent  non-rotating  part  the  components  in  the 
directions  of  any  instantaneous  position  of  these  axes  are  constant.  Let 
^epf) » > ®pf  y be  t*le  comPonents  of  the  local  propulsive  force  referred  to 

the  axes  fixed  (virtually)  in  the  rotating  part.  Then  our  assumptions  are 

constant  on  the  rotating  part  (7-91) 


constant  on  the  adjacent  non-rotating  part.  (7-92) 


The  propulsive  force  on  the  non-rotating  part,  referred  to  the  no-def ormation- 
body-fixed  axes  is  obtained  by  premultiplying  the  last  expression  by  n and  so 
will  not  have  constant  components.  We  are  able  to  avoid  greater  complication, 
such  as  that  of  Ref  4,  because  we  again  have  some  local  constraint  on  the  modes  - 
the  assumption  of  almost  rigidity  for  the  rotating  part  produced  the  conditions 
(7-4)  and  (7-10),  though  different  from  that  of  Ref  2;  while  in  Ref  4 there  was 
no  constraint  at  all  on  the  modal  matrix  function  R . 
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* Body-fixed  axes  are  axes  which  are  fixed  in  a small  material  portion  of  the 
body. 
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The  local  propulsive  force  vector  referred  to  the  no-deformation-body- 
fixed  axes  is  therefore  (of  equation  (7-8)) 
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where  the  matrices  P ^ are  analogous  to  the  B . of  equation  (7-8);  te  they  are 
symmetric  and  such  that 


>> 

pf 

ftp) 

pf 

,<P> 

l*PfJ 


(7-93) 


(11)  element  of  = ith  element  of  |e 


[■ 


,(“) 

Pi 


(p) 

f(p) 

(P)l 

pf 

fpf 

Vj 

(p) 

f (p) 

R(P)1 

pf 

pf 

8pfJ 

of 

e(n) 

P 

(see 

Hr 

J V 


a2n. 


p"  3x 


pf 


2 

rnt 


) t ,3y  c 
P pf  Pf 


etc 


(7-94) 


P = 


[•{? 

f(n) 

Pf 

nf 

an 

3x  , 
Pf 

- 

f(p) 

pf 

•S] 

[•£ 

f (n) 
Pf 

< 

anf  i 

aV 

[•S’ 

f(p) 

pf 

[•* 

,(n) 

Pf 

’’f 

anf 

3v 

f (p) 
pf 

3fL 


p t 3x  , 
P Pf 


P R 

J V 


an. 


_ 3y  , 
P pf 


°1R 

J Pnfc 


an. 


. 3z  , 
P Pf 


(7-95) 


analogous  to  the  expression  for  B (equation  (7-14));  then,  for  any  vector 

|f,  n C}, 
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and  so  (7-93)  can  be  rewritten  as 
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A further  condition  on  the  propulsive  force  components,  which  we  will 
apply,  is  that,  in  the  unperturbed  state,  the  total  propulsive  force  due  to  each 
rotating  part  will  have  constant  components  referred  to  the  no-deformation-body- 
fixed  axes,  and  so 

• t V'(p)rr„n 


*T  r 


(7-98) 


The  total  region  where  the  propulsive  force  acts,  denoted  by  (p+),  is  the 
rotating  part,  plus  a bit  more,  but  over  (p+)  - (p)  this  condition  is  satisfied 
as  a consequence  of  (7-92).  From  the  form  of  Rppt  ( of  equation  (2-4))  we 
then  have 
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(7-100) 


Finally,  we  assume  that  the  moment  of  the  propulsive  forces  on  the  rotating  part 
is  about  the  axis  of  rotation  and  so 
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Now  the  rotating  part  is  axisymmetric  and  so  any  section  of  it  normal  to 
the  axis  of  rotation  must  have  n-fold  symmetry,  about  the  point  of  intersection 
of  the  axis  and  the  section,  where  n is  an  integer  greater  than  or  equal  to  3. 
Consequently. 
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we  find  the  sums  of  the  left  hand  sides  of  these  expressions  are  both  zero,  and 
so  we  have 
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x any  function  of  x^p^ 
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These  relationships,  combined  with  (7-101),  also  enable  us  to  deduce  that 
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This  is  not  all  that  we  can  deduce  from  (7-102)  and  (7-103).  We  have 
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and  so  we  find  that 
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7.4.1  Generalised  propulsive  forces 

The  column  vector  of  the  generalised  propulsive  forces  is 
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Now,  of  equations  (7-16)  and  (7-97) 
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In  the  Appendix  we  established  a set  of  relationships  of  the  form 
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and  so  (of  equation  (7-95),  (7-14)  and  (7-18)) 
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where 
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etc. 
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,,  using  (7-110),  (7-111),  (7-107),  (7-105)  and  (7-106),  we  find  that 
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and  so 
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yl  , and  so,  for  example 
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Similarly, 
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etc.  In  this  last  expression,  the  row  matrices  are  the  transposed  columns  of 
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where  (c/  (7-15)) 
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Equations  (7-118)  and  (7-119)  could,  of  course,  as  a consequence  of  (7-96),  have 
been  written  in  terms  of  just  three  sums: 
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(p) 
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but  all  comeliness  is  then  lost.  We  will  not  transcribe  the  other  equations 
similar  to  (7-119)  but  merely  note  that 
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• (7-123) 

To  complete  the  exposition  of  the  second  terra  in  the  expression  (7-112),  for  the 
generalised  propulsive  forces,  we  require  also  the  following,  easily  obtained 
relationship: 
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For  the  first  term  in  (7-112)  we  require  the  following  relationship  which  is 
obtained  from  equation  (7-97)  using  (7-95),  (2-4),  (7-105),  (7-106)  and 
equations  (A-5)  to  (A-14) 
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(see  equations  (3-4)  and  (4-7)  for  the  definition  of  the  J and  D matrices). 
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The  other  term  in  (7-112)  involves  the  summation  of  the  local  propulsive 
forces,  acting  on  the  non-rotating  parts  of  the  aircraft,  multiplied  by  the  modal 
displacements.  If  we  write*  (cf  (7-92)) 

r-i(p+)-(p)  x t t rcnfl 


RW  , e(P 

'V  ? 


(7-126) 


where,  from  (7-92),  P*  is  constant,  then  we  require  (see  equations  (7-97), 

(7-95)  and  (7-93))  the  following  additional  relationships  which  are  easily  deduced: 
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* The  minus  sign  has  been  included  in  this  definition  for  consistency  with  our 
definition  of  the  propulsive  contributions  to  the  generalised  forces  as 


■i(--U*EpuV-f 


The  elements  of  P^  will  then  be  seen  to  be 


the  parts  of  the  Pjf  arising  from  the  propulsive  forces  on  the  non-rotating 
part  of  the  aircraft. 
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Consequently 
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Incidentally  we  note,  from  the  Appendix  (equations  (A-15)  to  (A-22))  and 
equation  (7-120)  that 
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where  |e^^  f^^E  ^pf^}  are  t*le  comPonents»  the  unperturbed  state,  of  the  local 
propulsive  force  referred  to  the  no-deformation-body-fixed  axes  ( of  equation 
(7-92)).  Also,  from  equations  (7-10)  and  (7-16),  on  the  rotating  part 
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and  so,  if  ve  write 
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(7-132) 


(7-133) 


then 
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(7-135) 


which  shows  that  C (R2“Rj)  is  akin  to  the  modal  slope  matrix  of  earlier 

work  (Ref  2,  equation  (11))  but  for  a point  on  the  rotating  part  rather  than  the 
origin  of  the  body-fixed  axes. 

Putting  all  the  results  together,  we  therefore  find  that 
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where  is  given  by  equation  (7-127),  and  as  a consequence  of  (7-92)  is  not 

time  dependent. 

7.4.2  Overall  propulsive  forces 
We  have 
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These  two  expressions  bear  a prima  facie  resemblance  to  the  expressions  for  the 
generalised  propulsive  forces  (equation  (7-112)).  However,  while  the  unit  matrix 
can  be  considered  as  a particular  case  of  R - it  satisfies  the  condition  (7-4) 
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cannot  since  the  x 
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etc  are  functions  of  the 


generalised  coordinates  q^  (equation  (7-17)).  We  can  therefore  write  down 
immediately,  by  comparison  with  (7-136),  the  overall  propulsive  forces  as 
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and  is  a constant  ( cf  (7-92)). 
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The  moments  about  the  origin  of  the  no-deformation-body-fixed  axes  require 
further  analysis.  Now,  on  the  rotating  part  y^  *8  8iyen  by 

equations  (7-17)  and  (7-18)  while  for  |epn^  fp1^  g^n^j  we  turn  to  equation 
(7-97).  By  comparison  with  (7-16)  it  is  clear  that 
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Thus  we  immediately  have  the  following  relationship  which  we  require 
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In  addition,  using  (7-88)  and  (7-99) 
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while,  from  (7-105)  to  (7-107),  (7-110)  and  (7-111), 
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With  a little  manipulation  one  can  then  show  that 
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it  is  easily  shown,  in  a similar  manner  to  that  used  for  the  contributions  to  the 
generalised  forces  from  the  propulsive  forces  on  the  non-rotating  parts,  that  the 
same  forces  produce  moments: 


(7-151 ) 

The  expressions  for  the  total  moments  are  therefore  given  by  adding  this  equation 
to  (7-148).  Because  of  their  length  we  will  not  give  them  here. 

8 STRUCTURAL  FORCES 

Resulting  from  such  things  as  the  stresses  in  the  material,  friction 
between  different  components  of  the  structure,  etc,  there  will  be  a force  on  any 
particle  of  the  aircraft  which  can  be  called  a structural  force.  Calling  the 
components  of  this  structural  force  vector  (e^n\  f^n\  g , referred  to  the 

VS  S S ! 
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This  assumes  a linear  relationship  between  stress  and  strain,  is  consequent  upon 
the  linear  relationship  between  a particle's  cartesian  coordinates  and  the 
generalised  coordinates  (equation  (2-1)),  and  approximates  to  any  contribution 
not  resulting  from  a stress  by  a contribution  of  the  above  form.  The 
coefficients,  eg^  etc,  may  be  differential  operators. 

The  structure  cannot  exert  any  overall  force  or  moment  on  itself,  and  so 
it  is  easily  shown  (of  Appendix  A of  Ref  2)  that  the  following  conditions  must 
be  satisfied: 


The  last  of  these  equations,  (8-6),  comes  from  equating  to  zero  the  second  order 
terms,  in  the  generalised  coordinates,  in  the  expression  for  the  moment  about, 
say,  the  origin  of  the  no-deformation-body-fixed  axes.  If  we  had  included  some 
second  order  terms  in  the  expression  for  the  structural  force  (equation  (8-1)) 
this  condition  (8-6)  would  no  longer  hold  though  the  other  four  would. 
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Generalised  structural  forces 


The  column  vector  of  the  generalised  structural  forces  is 


4<»> 


• E- 


‘Tfesfl  + rTFs1  •••  esjfqi 


(8-7) 


8.2  Overall  structural  forces 

The  overall  translational  forces  and  moments  will,  as  already  stated,  be 
zero  whatever  the  perturbation  and  datum  state.  This  fact  has  indeed  already 
been  used  to  derive  the  conditions  (equations  (8-2)  to  (8-6))  which  are 
satisfied  by  the  coefficients  in  the  expression  (8-1),  for  the  structural  force 


vector. 


EQUATIONS  OF  EQUILIBRIUM 


The  equations  of  motion  will  of  course  be  satisfied  whatever  the  perturba- 
tion from  the  datum  state  and  in  particular  when  the  perturbations  are  zero. 

The  equations  then  obtained,  we  describe  as  the  equations  of  equilibrium  or 
datum  motion  equations.  In  this  state  the  effective  forces  are  all  zero 
( of  equations  (4-16)  to  (4-18),  (7-71),  (7-86)  and  (7-90)),  and  so  the  equations 
merely  state  that  the  sums  of  all  the  contributions  to  the  generalised  and 
overall  forces  are  all  zero.  Thus  we  have 
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(9-1) 


where,  from  (5-11),  (5-12)  and  (5-13) 
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| and  from  (8-7) 

S 

i 


- [(Ei>£]  ■ 


I' 


sf 

f8f 

8sf 


(9-5) 


The  contribution  to  the  propulsive  forces  from  the  non-rotating  part,  denoted 
above  (equation  (9-4))  by  all-embracing  symbols,  is,  from  (7-126),  (7-140)  and 
(7-149) 
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where  the  local  propulsive  forces  satisfy  (7-92). 
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9. 1 Other  equilibrium  states 

Once  non-linear  terms  are  introduced  into  the  equations  of  motion  there  is 
always  the  possibility  of  equilibrium  states,  other  than  the  datum  state  when 
all  the  perturbations  are  zero.  Any  such  equilibrium  states  will  be  given  by 
solutions  of  the  equations  of  motion  which  satisfy  the  condition  that  all  the 
generalised  coordinates  are  constant.  Thus,  for  example,  if  we  had  no  deforma- 
tional  freedoms  (a  rigid  aircraft),  possible  equilibrium  states  are  given 
approximately  by  the  solutions  of  the  equation  ( of  equations  (5-14),  (5-15), 
(6-9)  and  (6-10)): 
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(9-7) 

where  the  subscript  qs  denotes  the  quasi-steady  values.  The  above  equation 
neglects  any  third  or  higher  order  terms  in  the  generalised  coordinates.  No 
doubt  this  often  will  not  be  good  enough  for  the  determination  of  any  equilibrium 
state  which  is  not  the  datum  state.  In  such  circumstances  a better  approxima- 
tion to  the  local  aerodynamic  force  vector  than  (5-8)  will  pretty  certainly  be 
required;  and  so  an  entirely  different  approach,  such  as  solving  the  aerodynamic- 
gravitational-propulsive-structural-dynamic  problem  as  one  problem  rather  than 
doing  aerodynamic  calculations  etc  first  to  provide  data  for  the  dynamic  problem. 

It  should  also  be  noted  that  if  an  autonomous  system  such  as  ours  has  two 
equilibrium  states,  one  of  which  is  stable  then  the  other  one  will  be  unstable. 
However  one  can  have  two  equilibrium  states  both  of  which  are  unstable*. 

* At  an  equilibrium  point  the  eigenvalue  equation  has  the  form  X^n+ajX^n  1 + ... 
+ a2n * 0 • there  are  just  two  such  points  then  a2n  at  one  will  have  the 
opposite  sign  to  a2n  at  the  other.  Therefore  all  the  eigenvalues  at  both 
points  cannot  all  have  negative  or  zero  real  parts. 
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10  PERTURBATION  MOTION  EQUATIONS 

The  equations  of  motion  for  perturbations  from  the  datum  state  clearly 
follow  from  (4-1),  (4-17)  and  (4-18)  when  all  the  contributions,  to  the  various 
constituents  of  the  generalised  and  overall  forces,  which  are  independent  of  the 
generalised  coordinates  are  omitted  (af  (9-1)).  In  the  previous  sections  we  have 
obtained  expressions  for  the  various  elements  in  these  equations.  The  relevant 
contributions  to  the  effective  forces  are  displayed  in  equations  (4-16)  to  (4-18), 
(7-71),  (7-86)  and  (7-90);  while  those  to  the  applied  forces  (generalised  and 
overall)  are  found  in  equations  (5-11),  (5-14),  (5-15),  (6-8)  to  (6-10),  (7-136), 
(7-139),  (7-148),  (7-151)  and  (8-7).  There  is  of  course  no  structural  contribu- 
tion to  the  overall  forces.  Little  would  be  achieved,  other  than  indigestion, 
if  we  were  to  gather  together  here  all  these  various  contributions,  except  in  the 
case  of  small  perturbations  which  we  will  consider  in  a moment. 

However  it  is  perhaps  instructive  at  this  stage  to  consider  briefly  the 
sort  of  behaviour  that  might  be  predicted  by  our  non-linear  equations.  For 
example  take  a one  degree  of  freedom  system  whose  equation  is 

••  • 2 

q + bq  + cq  + dq  =0  b,  c > 0.  (10-1) 

It  has  two  equilibrium  states: 

(i)  the  origin  q = 0 which  is  a point  of  asymptotic  stability;  and 

(ii)  the  point  q = -c/d  which  is  a saddle  point  and  therefore  unstable. 

The  trajectories,  in  the  (q,q)  plane  then  take  the  form  shown  in  Fig  1 (see 
Ref  13,  p.62).  It  will  be  seen  that  the  region  of  asymptotic  stability  is 
bounded . 

A second  example  is  the  one  degree  of  freedom  system  whose  equation  is 

q + hq^  + cq  - 0 c > 0.  (10-2) 

In  this  case  we  have  only  one  equilibrium  state  (q  - 0)  but  the  region  of 
stability  is  still  bounded  {af  Ref  14,  p.104)  - see  Fig  2.  The  same  sort  of 
thing  is  true  for  the  equation 


q + kqq  + cq  - 0 


c > 0.  (10-3) 
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(see  Fig  3).  In  both  these  cases  ((10-2)  and  (10-3))  the  datum  (equilibrium) 
state  is  stable  but  not  asymptotically  stable,  but  it  is  possible  to  have  a 
bounded  region  of  stability  even  when  there  is  only  one  equilibrium  stat^  and 
that  is  asymptotically  stable.  This  is  illustrated  by  the  example 

q + (bq  + cq)^l  + | qj  = 0 k,b,c  > 0,  (10-4) 

for  which  the  region  of  asymptotic  stability  is  bounded  by  q = - c/k  (see  Fig  4). 

Another  possibility,  often  encounted  in  physical  systems,  is  the  occurrence 
of  limit-cycles.  Linear  theory  then  predicts  the  equilibrium  state  to  be 
unstable  but  non-linear  theory  shows  that  it  is  stable  in  that  any  path  starting 
in  the  vicinity  of  the  equilibrium  point  tends  ultimately  to  a closed  curve 
surrounding  that  point.  This  path  is  then  a stable  limit-cycle*.  Such  solutions, 
in  the  one  degree  of  freedom  case,  will  not  be  predicted  by  our  approximation. 

It  is  necessary  to  include  third  order  terms  at  least.  With  more  than  one 
degree  of  freedom  one  cannot  say  that  no  limit  cycles  are  possible  with  present 
approximation  but  here  again  one  would  intuitively  think  that  a better  approxima- 
tion is  required. 

To  sum  up  therefore,  one  may  say  that  the  non-linear  equations,  as 
developed  here,  may  be  useful  to  indicate  whether  there  is  a boundary  to  the 
region  of  stability,  but  probably  one  of  little  value  in  determining  whether  an 
unstable  equilibrium  state  is  'really'  unstable  or  possesses  a stable 
1 imit-cycle . 

10.1  The  linear  approximation 

When  we  collect  together  from  sections  4.1,  4.2,  5.1,  5.2,  6.1,  6.2,  7.2, 
7.3,  7.4.1  and  7.4.2  the  various  contributions  to  the  applied  and  effective 
generalised  and  overall  forces  we  find  that  for  small  perturbations  the  following 
equation  must  be  satisfied: 


{[A..]D2  . [J..]D  • [V..]  . [G,.]  . [P;j]  * [EU]  - [Qij]} 


ln+6 


(10-5) 


* There  is  also  the  possibility  of  unstable  limit-cycles.  The  equilibrium  point 
is  then  asymptotically  stable  - any  path  starting  within  the  limit-cycle  will 
finish  up  at  the  equilibrium  point.  Semi-stable  limit-cycles,  stable  on  one 
side,  unstable  on  the  other,  should  not  occur  in  practice. 


i 

t 

j 


■j 

] 


(10-7) 


(iii)  Ponderous  stiffnesses 
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(10-12) 


These  equations  differ  somewhat  from  what  might  be  considered  to  be  the 
corresponding  case  among  those  considered  in  Ref  2 (see  Ref  2,  Table  6 "Body- 
fixed  axes,  free-free  modes,  displacement  body  freedoms").  The  difference  is 
not  solely  due  to  the  different  propulsive  force  model  (section  7.4)  and  to  the 
inclusion  in  the  present  development  of  the  effect  of  a rotating  engine  which 
contributes  the  ponderous  damping  and  ponderous  stiffness  terms  (equations  (10-6) 
and  (10-7)).  As  mentioned  in  the  Introduction,  the  basic  cause  is  a slightly 
different  representation  of  the  deformations.  Thus  from  equations  (12)  and 
(103)  to  (107)  of  Ref  2 it  is  seen  that  in  that  derivation 
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This  agrees  with  the  present  representation,  equation  (2-2),  as  far  as  the  first 

order  terms,  and  the  agreement  is  complete  if  Rq  and  are  zero.  These  two 

matrices  - matrices  of  modal  displacements  and  slopes  at  the  reference  point  - 

appeared  in  Ref  2 because  of  the  indirect  derivation  from  the  equations  of  motion 

for  the  encastr£  modes  case.  If  we  put  R.  and  P zero  in  Table  6 of  Ref  2, 

r 0 q 

the  gravitational,  structural  and  aerodynamic  matrices  agree*  with  those  given 
above  (equations  (10-9),  (10-11)  and  (10-12)).  Furthermore  if  we  put  R£  = Rj 
and  assume  all  the  propulsive  force  is  on  a non-rotating  part  - this  makes  the 
two  propulsive  force  models  the  same  - then  the  only  non-zero  term  in  (10-10)  is 
the  (-  R)  and  there  is  agreement  with  the  propulsive  matrix  of  Ref  2,  Table  6, 
when  R„  = P = 0 . 


* The  ponderous  inertia  matrices  agree  in  any  case. 
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the  derivatives  of  n,  nf 


n = n(xp,  yp,  Zp) 


n. 


where,  cf  equation  (7-7),  with 


n(xpf’  ypf*  ZpP 


2 2 2 
x + y + z 


2 2 . 2 
s * x + y 


II  (x 


,y,z)  and  its  required  derivatives  are  given  by 


a_n 

3x 


an 

ay 
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(A-0 
(A- 2) 


(A-3) 

(A-4) 


n = 

" x/r 

- y/s 

- xz/rs 

— 

y/r  z/r~ 

x/s  0 

-yz/rs  s/r 

J 

f (r2-x2)/r3 

-*vA3 

/ 3 “ 

- xz/r 

xy/s3 

,2A3 

0 

[Z(x2s2-y2r2)/(rs)3 

xyz (r2  + s2)/ (rs) 

3 xz2/r3s 

-xy/r3 

Cr2-y2)/r3 

-yz/r3 

2/  3 
- x / s 

-xy/s3 

0 

xyz(r2  + s2)/ (rs) 

z(y2s2  - x2r2)/ (rs) 

3 2/3 

yz  /r  s 

(A-5) 


(A-6) 


(A-7) 
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Appendix 


an 

3z 

r / 3 

- xz/r 

-yz/r3 

2/3 
s /r 

(A-8) 

0 

0 

0 

• 

2/  3 

- xs  /r  s 

-ys2/r3s 

/ 3 

- zs/r 

t 

$ ■ 


3x(x2-r2)/r5 

y(.2-3x2)/*5 


y(3x2-r2)/r5 

-3xy2/.5 


x(3x2  -r2)/i 
0 


xz 

(r.)1 


f3r2.2(r2*.2)  J -JT*  (AV*f)  2 2 4 l -HtrV-«2(t2*3.2)l 

1 - x2(3r4*2rV*3.4)J  0^7  1 - x2(3r4  ♦ 2rV  ♦ 3.4)  1 r. 


(A-9) 


a2n 

7 


3y(y2-r2)/r5 
x(3y' 


*C3y2-r2)/rS 
0 


2-.2)/.5 


XX 


x(3y2  -r2)/t5 
3x2y/.5 

frV(r2*.2)  1 .,J«  [3r2,y  **2)  .1  -f-y  |r2»2  -y2(r2  + 

{ - y2(3r4  + 2r2.2  + 3.4)J  (r«)5  [ - y2(3r4  ♦ 2rV  ♦ *4>j  r» 


3»2)} 


32n 

3z2 


x(3z 


2 - r2)/r5  y(3z2  - r2)/r: 


L 3sxz/r 


3syz/r 


-3zs2/r5  • 

0 

s(3z2  - r2)/r5_l 


(A-10) 


(A— 1 1 ) 


a2ti 

3x3y 


y(3x2-r2)/r5 


-*(3y 


2 »2)/.5 


x(3y2-r2)/r5 

- y(3y2  - 2.2)/*5 


yt  JrV  I **  Jr2*2(r2*«2>  l 

(r«)S  j - x2(3r4*2rV  + 3»4)j  <r»)5  1 - y2(3r4  ♦ 2r2*2  + 3*  )J 


3xyr/r5 

0 


-2«-2V)/rV 


- xyr  (3*  +r 


J 

(A-l  2) 
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a2n 

3y3z 


3xyz 

0 


A5 


W-rVr; 

0 


xy(3s2  - r2)/r3s  - (s/r3)  + y2  (3s2  - r2)/r5s  yz(3s2-r2)/r 


y(3s2  -2r2)/r5 
0 

2 2X  /_5 


(A- 13) 


and 


32n 

3z3x 


z(3x2  -r2)/r! 
0 


3xyz/:‘ 

0 


- x(3s2  - 2r2)/r5 


- (s/r3)  + x2(3s2  - r2)/r5s  xy (3s2  - r2)/r! 


s xz 


(3s2  -r2)/r5s 
(A— 14) 


t an  t 

Since  IUI  = I , the  products  -z—  II  etc  are  skew  symetric  and  we  find 

dX 


that 


^nT  - 

3x  a 


(A- 15) 


where 


|inT 

3y 


= Al 


|£nT 

3z 


= A. 


a 


JJL 

2 

rs 


b 


2 

r s 


(A-l 6) 

(A— 1 7) 

(A-l 8) 

(A-l 9) 

(A-20) 
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Append ix 


Consequently 


with 


fin 


T 


p 


[a  b 


- - 

yz/rs2 

- zx/rs2 

0 

• 

X 

- zx/r2s 

- yz/r2s 

CM 

U 

to 

• 

y 

y/rs 

- x/rs 

°j 

• 

Z 

(A-2 1 ) 


(A-22) 


a 


It  is  easily  seen  also  that 


_an  3nT 

9x  3x 


an  an^ 

3x  3y 

an  an^ 

ax  az 


an  an^ 

ay  3y 


(A-23) 


-A2 

a 

(A-24) 

‘ Vfc 

(A-25) 

- A A 

(A-26) 

a a 

(A-2 7) 

The  transpose  of  the  square  matrix  that  appears  in  (A-22),  ie  the  matrix 

T T 

j a b a (,  has  been  denoted  by  C in  the  main  part  of  the  report  for  the 
case  when  the  arguments  are  x ^ etc  (of  equation  (7-15)). 

A further  relationship,  which  can  be  easily  verified,  which  we  will  use  is 


|A  r 


1 

0 

LoJ 


A 
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3y  3z 


JT  " - V. 


3n  3lr 

3z  3z 


- a: 


and  so 


il  nT 
2 
3x 


l!lLnT 

3x3y 


3H  nT 


3A 

< 

9x 


a + A2 


3x3z 


i^nT 

3y" 


£lLnT 

3y3z 


i!inT 

a 2 

3z 


2k 

— — + A A, 
3y  a b 


2 k 

•—2-  + A A 
3z  a c 


3A,  •> 

V1  + AJ 

3y  b 


3 A., 

+ A,  A 
3x  b a 


3AC 

— — — + A A 
3x  e a 


3Aj 

3z 


+ Vfl 


3A 

—2-  + A A, 
3y  c b 


3AC  2 

— + AZ 
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(A-28) 

(A-29) 

(A-30) 

(A-31) 

(A-32) 

(A-33) 

(A-34) 

(A-35) 
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A,  etc 
9 


E. 

1 


E.  . 
ij 

E(y),E(pv) 
r * r 


I 

n 


J. 
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LIST  OF  SYMBOLS 

skew-symmetric  matrices  involving  etc  (see 

equation  (2-9)) 

ponderous  inertia  coefficient 

see  equation  (2-10) 

see  equation  (4-9) 

differential  operator  d/dt 
- diag{$  6 

generalised  structural  force 
structural  stiffness  coefficient 

coefficient  matrices  in  second  order  expression  for  aero- 
dynamic forces  (see  equations  (5-3)  and  (5-7)) 

generalised  gravitational  force 

gyrostatic  force  (equation  (4-4)) 

gravitational  stiffness  coefficient 

unit  matrix  (generally  3x3) 

2 

* - 6mA  (matrix  of  moments  and  products  of  inertia) 

Xf 

a certain  coupling  force  between  the  rotational  body  freedoms 
and  the  deformational  freedoms  in  Lagrange's  equations 
referred  to  a non-inertial  frame  (see  equation  (4-5)) 

a certain  matrix  formed  from  the  elements  of  {$  0 i|»)  - see 
equation  (3-4) 

ponderous  damping  coefficient 


K, 


L 


L. 

l 


L-.L-  etc 


M*,M*  etc 
x 9 

N 


coefficient  matrices  in  second  order  expression  for  aero- 
dynamic forces  (equation  (5-7)) 

rolling  moment 

rolling  moment  coefficient  (aerodynamic) 
aerodynamic  rolling  moment  coefficients 
pitching  moment 

pitching  moment  coefficient  (aerodynamic) 
aerodynamic  pitching  moment  coefficients 
yawing  moment 
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LIST  OF  SYMBOLS  (continued) 


N. 

1 


yawing  moment  coefficient  (aerodynamic) 


' ) factor'  in  axes  transformation  matrix  (see  equation  (2-5)) 


P.  . 
ij 


generalised  propulsive  force 

matrix  of  contributions  to  generalised  propulsive  force 
resulting  from  the  propulsive  forces  on  the  non-rotating  parts 
of  the  aircraft 

propulsive  stiffness  coefficient 


generalised  force 


R 


matrix  relating  angular  velocities  and  orientation  (see 
equation  (3-3)) 

modal  matrix  (see  equation  (2-1)) 


R1  ,R2 


values  of  R at  two  points  on  axis  of  symmetry  of  a rotating 
part 

'tj>  factor*  in  axes  transformation  matrix  (equation  f 2—4 ) ) 


S,  Scf> 

f 


axes  transformation  matrix.  Absence  of  a subscript  means 
the  arguments  are  9,  (see  equations  (2-3)  and  (2-8)) 


V 


0 


centrifugal  potential  function  (see  equation  (4-2)) 


V.  . 

1J 


w 


X,Y,Z 
X . , Y . ,Z  . 

l l l 


XA’X* 

YA’Y* 

zx  z*j 


etc 


ponderous  stiffness  coefficient 

kinetic  energy  relative  to  the  frame  of  reference  (see 
equation  (4-3)) 

overall  force  resolutes 

overall  force  resolute  coefficients  (aerodynamic) 


aerodynamic  force  resolute  coefficients  - they  may  be 
differential  operators 

*ij>  factor*  in  axes  transformation  matrix  (equation  (2-6)) 


a . ,b . ,c . 
ill 

elements  of 

R (equation  (7-134)) 

,g 

e. ,e  ,e  " 
i x 4> 

f . ,f  ,r 

components 

of  local  force  vector 

.etc 

local  force 

vector  coefficients  - they  may  be  differential 

1*  x*  $ 

8i,8x*S. 

operators 
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fi*f* 

8x,8i 


s 

m 

6m 

n 

p.q.r 


etc 


p("0>,  ,«». 


.(pO) 


"1 

2 

r 

2 

s 

t 

u,v,v 


u ,v  ,w 
m’  m m 


x,y,z 
x 1 * y 1 * z 1 


x ,y  ,z 
P P P 


Bf,\’BK2’BK3 

C 

P 

Prp2.P3 


LIST  OF  SYMBOLS  (continued) 


local  aerodynamic  force  vector  coefficients  - they  may  be 
differential  operators 


acceleration  due  to  gravity 

first  column  of  S^,  (see  equation  (6-5)) 
f 

second  column  of  S^,  (see  equation  (6-6)) 

f 

third  column  of  S^  (see  equation  (6-3)) 
f 

mass  of  aircraft 
mass  of  a particle 

{number  of  def ormational  degrees  of  freedom. 

Order  of  symmetry  of  section  of  rotating  part 

angular  velocity  resolutes 

angular  velocity  of  rotating  part 

see  equation  (7-65) 

generalised  coordinate 

2 2 2 

= x + y + z 

2 2 

* x + y 

time 

linear  velocity  resolutes 
particle  velocity  resolutes 

particle  position  resolutes 

position  resolutes  for  the  origin  of  the  no-deformation-body- 
fixed  axes  relative  to  the  origin  of  the  constant-velocity 
axes 

see  equations  (7-1)  and  (7-2) 


matrix  relating  modal  matrix  for  a general  point  on  the 
rotating  part  with  that  for  a point  on  the  axis  of  symmetry 
(see  equations  (7-10),  (7-12),  (7-14)  and  (7-16)) 

see  equation  (7-8)  et  eeq 

see  equation  (7-15) 
see  equation  (7-95) 
see  equation  (7-94) 

see  equation  (7-150) 
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LIST  OF  SYMBOLS 


(continued) 


a,b,a 

see 

aPfV,e 

pf 

the 

are 

aa  ab 

aa  ] 

ba  bb 

ba 

the 

aa  ab 

aa  \ 

equations  (A-18)  to  (A-20) 

T 

columns  of  C , which  are  a,  b,  a when  the  arguments 
x , etc 
Pf 

T 

elements  of  BP  (see  equations  (7-115)  and  (7-116)) 


e_ 


n 


*f 

% 

thus 


n 

a.e.Y 

£,n,c 


diag  (c,n,c) 


E 

E 

I 

E 


(p) 

(p+) 


(p+)-(p) 


angle  of  inclination  in  datum  motion 

product  of  *9  and  $ factors'  in  axes  transforn ation  matrix 
for  transformation  from  the  no-deformation-body-f ixed  axes 
to  the  rotating  part-'fixed'  axes  (see  equation  (7-7)) 

angle  of  bank  in  datum  motion 
nose-azimuth  angle  in  datum  motion 

are  the  orientation  angles  of  the  constant-velocity  axes 
relative  to  the  normal  earth-fixed  axes 

see  equation  (A-5) 

see  equations  (7-50)  and  (7-50b) 

elements  of  arbitrary  vector 

orientation  angles  of  no-deformation-body-fixed  axes 
relative  to  constant-velocity  axes 

the  diagonal  matrix 


indicates  summation  for  all  the  particles  of  the  aircraft 

indicates  summation  over  the  rotating  part 

indicates  summation  over  the  total  region  where  the  propulsive 
force  acts 

(p+) 


COO 
0 n 0 

o o d 


96 


LIST  OP  SYMBOLS  (concluded) 

DRESSINGS 

(i)  Subscripts  i 

A quantity  is  relative  to  or  about  the  origin  of  a particular  set  of  axes  j 

r 

where 

c denotes  the  constant-velocity  axes 

n denotes  the  no-deformation-body-fixed  axes 

p denotes  axes  'fixed'  in  the  rotating  part  j 

f indicates  the  value  during  the  datum  motion 

g indicates  gravitational  i 

p indicates  propulsive  j 

j 

s indicates  structural  j 

Occasionally  a subscript  is  added  to  an  array  to  show  that  all  the  elements  j 

of  the  array  should  have  that  subscript.  | 

! 

(ii)  Superscripts 

Bracketed  superscripts  denote  the  axes  of  resolution  with  the  same 
significance  as  in  (i)  above.  Also 

T denotes  the  transpose  of  a matrix 

+ indicates  a contribution  from  the  non-rotating  parts  of  the  aircraft  where 

propulsive  forces  act 

(pO)  see  equation  (7-65) 

(iii)  Suprascripts 

• (dot)  denotes  derivative  with  respect  to  time 

- (bar)  denotes  typical  or  total 
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